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We develop a theory of extrinsic spin currents in semiconductors, resulting from spin-orbit coupling 
at charged scatterers, which leads to skew scattering and side jump contributions to the spin Hall 
conductance. Applying the theory to bulk n-GaAs, without any free parameters, we find spin 
currents that are in reasonable agreement with recent experiments by Kato et al. [Science 306, 
1910 (2004)]. 



Generating and manipulating non-equilibrium spin 
magnetization by electric fields is one of the most desir- 
able goals of semiconductor spintronics, because electric 
fields have potentialities for accessing individual spins at 
nanometer scale. Spin-orbit (SO) coupling is a mecha- 
nism for achieving this goal. It has the prominent conse- 
quence of the spin-Hall effect (SHE), where an electric- 
current can induce a transverse spin current and a non- 
equilibrium spin accumulation near sample boundaries. 
Recent observations of this effect are important experi- 
mental achievements 0, . Theoretically, two different 
mechanisms of SHE were proposed. The extrinsic mech- 
anism 0, 0, is based on spin-dependent scattering of 
electrons by impurities and is mostly due to Mott skew 
scattering |Q| . An intrinsic mechanism has also been pro- 
posed, which is related to the concept of dissipationless 
spin currents in a perfect crystal 0] • 

The theory of spin transport in media with SO cou- 
pling is rather intricate and includes all problems inher- 
ent in the theory of anomalous Hall effect ( AHE) , which 
has a long history; for reviews see 0,0]. Also, the pre- 
cise definition of spin currents is still under dispute be- 
cause of the problems stemming from spin nonconserva- 
tion in media with SO coupling. If spin-relaxation rates 
are small, a spin current with non-zero divergence can 
lead to spin accumulations, which are experimentally ob- 
servable quantities 0, 0, ^| . However, because spin- 
currents are even with respect to time inversion, they 
do not necessarily vanish in thermodynamic equilibrium, 
hence their relation to spin transport and spin accumu- 
lation is far from obvious 

These problems inherent in the spin-transport theory 
make identification of physical mechanisms underlying 
the SHE observed in Refs. 0] and 0rather challenging. 
On the one hand, Wunderlich et al. 0] observed a strong 
SHE in two-dimensional (2D) layers of p-GaAs and as- 
cribed it to the intrinsic effect because of the large magni- 
tude of the effect and large splitting of the energy spec- 
trum typical of heavy holes. On the other hand, Kato 
et al. 1] attribute their measurement of SHE in three- 
dimensional n-GaAs layers (2/xm thick) to the extrinsic 
mechanism. We believe that this is indeed the case as we 
explain in this work. Although an intrinsic spin Hall ef- 
fect, driven by the fc 3 -Dresselhaus SO coupling [l5|, could 
give rise to spin accumulation in this system, as proposed 
in Ref. Iltil its estimated size, when impurity scattering is 
taken into account, is an order of magnitude smaller than 




Figure 1: Spin-dependent scattering at an attractive impurity. 
We show the classical trajectories (solid lines), for a screened 
Coulomb potential and for strongly exaggerated spin-orbit 
coupling with A > and with quantization axis z. The skew- 
scattering current results from different scattering angles for 
spin- 1 and spin-| and leads to a positive spin Hall conductiv- 
ity, 0-$$ = —jss,y/Ex > 0. Further, we show the horizontal 
displacement due to the side jump effect (dashed lines) , con- 
tributing to the spin current with opposite sign. 

the observations. Further, because of the large sample 
thickness, a specifically 2D mechanism of spin accumu- 
lation proposed recently be Usaj and Balseiro 01 an d 
Nikolic et al. 0] and relying on the properties of near- 
edge states cannot play a role in the geometry of Ref. 1. 

In the following, we develop a theory of extrinsic spin 
currents in a 3D electron system. It results from intrin- 
sic SO coupling in the bulk crystal that produces a SO 
contribution to the impurity potential. (This effect can 
occur even when the host crystal is inversion symmetric, 
so that spin is conserved outside the radii of the impu- 
rity potentials.) We find that impurity scattering and 
SO interaction in n-GaAs are strong enough to support 
spin-currents that are in reasonably good agreement with 
findings by Kato et al. 0] without using any adjustable 
parameters. Although the sign of our result is opposite 
to that mentioned in Ref. [ij, there remains some uncer- 
tainty about the absolute sign in the experiments |l9j| . 
Since our result contains two contributions of opposite 
sign, and we have employed several approximations, we 
would believe the extrinsic mechanism to be relevant even 
if it turns out that the difference in the sign of the effect 
persists. 

We consider an electron Hamiltonian of form 

ff=il+y(r) + A<T.(kxVy), (1) 

2m* 

where V(r) is an electron potential energy that varies 
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slowly on the scale of the host lattice constant. In vac- 
uum, the last term of Eq. results from relativistic cor- 
rections in the Pauli equation and is known as Thomas 
term, with A = — h 2 /Arrive 2 ~ —3.7 x 10~ 6 A 2 , vacuum 
electron mass too, and velocity of light c. In direct gap 
cubic semiconductors like GaAs, a SO interaction of the 
same form develops in the framework of the k • p model 
due to the coupling of a s-type conductance band to p- 
type valence bands. This coupling is well established and, 
for asymmetric 2D systems, it leads to the Rashba term. 
For conduction band electrons in the 8x8 Kane model, 
one finds A = (P 2 /3) [1/JSg -l/(E + A f] in third- 
order perturbation theory, with gap Eq, SO splitting Ao 
between the J = § and J —\ hole bands, and a properly 
normalized interband matrix element P of the momen- 
tum |20j. All these parameters are large (i.e., of atomic 
scale) since they result from the strong crystal poten- 
tial. For example, for GaAs one finds A = 5.3 A 2 . Thus, 
the SO coupling in n-GaAs is by six orders of magni- 
tude stronger than in vacuum and has the opposite sign. 
This enhancement of SO coupling is critical for develop- 
ing large extrinsic spin currents. 

In the 3D problem that we consider, SO coupling comes 
only through the potential energy V(r) of Eq. |QJ, that 
can be either electron energy in an impurity center or its 
energy — e(E • r) in a driving electric field, where e is the 
charge of an electron, i.e., e < 0. We ignore here the 
SO /c 3 (Dresselhaus) correction to the electron energy in 
pure GaAs, which is small and is absent in the 8x8 Kane 
model. 

We now analyze the effect of the SO coupling in Eq. (JTJ 
on the scattering at impurities, which leads to the extrin- 
sic spin Hall effect. This comprises two contributions, 
one resulting from the skew-scattering at the impurities 
and the other from the shift of the scattered wave 
packet pi) (known as side jump contribution). These 
contributions are shown schematically in Fig. 1. 

The skew scattering contribution is calculated in the 
lowest order both in SO interaction and in l/k-pi, with 
Fermi momentum kp and mean free path £. To this end, 
we describe the system by a spin-dependent Boltzmann 
equation and a distribution function written as a 2 x 2 
spin matrix / = [/o(k) + </>(k)] 1 + f(k) • er, with equi- 
librium distribution function fg. In the following, we 
suppress the identity matrix 1. The collision integral on 
the right hand side of the Boltzmann equation has the 
form 



impurity potential, we can then write 



coll 



k': k'= 



hk da r 
to* d£l 



/>)"/>" 



where rii is the impurity density. The scattering cross 
section da/dfl is spin-dependent and mixes spin compo- 
nents of the incoming flux. In the general case, this spin 
dependence is rather complex [22j. However, it simplifies 
essentially when we expand Eq. J2J) in SO coupling and 
neglect the cross terms containing spin-dependent contri- 
butions of both da / dVt and /. For a central symmetric 



da r 
dh 
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where z9 = i^k' is the angle between k' and k, and 
n = k' x k/|k' x k| is the unit vector normal to the scat- 
tering plane. The coefficient /(i?) is the spin-independent 
part of the scattering cross section, while 5(i?) is the so- 
called Sherman function 0, 12^, E3 , which measures the 
polarization of outgoing particles scattered into direction 
k from an unpolarized incoming beam of momentum k'. 

To lowest order in the electric field, the left hand side 
of the Boltzmann equation equals eE ■ = (eh/m*)(E ■ 

k) where the isotropy of /o(k) was used. The Boltz- 
mann equation may then be solved with the following 
ansatz. First, the usual spin-independent term is set to 
0(k) = k • E Cfc [typically, C k = - (ehr/m*) ^ with 
transport lifetime r] . Second, we use f (k) = (E x k) D k 
for the spin polarization due to SO interaction. This 
structure is motivated by the physics of Mott scattering, 
namely, that the spin polarization is perpendicular to the 
scattering plane defined by incoming electrons that drift 
in the direction of — E and are scattered into k. Here, 
Cfc and D k are spherically symmetric functions of k. 

With this ansatz, we can evaluate the collision inte- 
grals on the right hand side of the Boltzmann equation. 
Integrating over the direction of k', dil (k') = dip d-d sin 
where ip is the the azimuth of k' in the plane perpendic- 
ular to k, and suppressing an overall factor of riihk/m* 
in Eq. we obtain (see Add. \K$ 

J dO(k') / (tf) (k - k') • (E Cfc + er x E D k ) 

= [k • E Cfc + k • (er x E)D k ] J dQ I (tf) (1 - cos tf) , 

(4) 

where the integral on the right is proportional to the 
inverse transport time r , and 

J dft(k') / (tf) S (i?) (er • n) (k • E - k' • E) C k 

= -i k • (cr x E) C k JdQ /(-!?) S{$) sin >&. (5) 

Since the left hand side of the Boltzmann equation 
only depends on the component of k along the electri- 
cal field, this must also be the case on the right hand 
side. Thus, the second term in Eq. 10} must cancel with 
Eq. determining D k — \^ k C k . We defined the trans- 
port skewness 
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which describes the effect of skew-scattering on the dis- 
tribution function and depends on the structure of the 
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scattering center and on the energy of the scattered par- 
ticle. Therefore, our ansatz is self-consistent, and the 
solution of the Boltzmann equation is 



/(k) = /o(*0 



E- 
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Cfc, 



(7) 



i.e., the components of f(k) are / M (k) 



(Exk) C t 



Now we calculate the contribution of skew-scattering 
to the spin current, jg S = n {a^ vo), with density n and 
with vo = ftk/m* (SO contributions to the velocity are 
analyzed below). We obtain 
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^3^/(k) = ^ £ ^(Jo),, (8) 
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where J = 2e J d 3 k (2ir)~ 3 (hk/m*) k-E Ck is the charge 
current in the absence of SO coupling, and summation 
over v is implied. Assuming low temperatures, we need 
to evaluate 7^ only near the Fermi energy Ep and we de- 
fined 7 = 7fe F [24| . If there are different species of impu- 
rity potentials, the weighted average of the corresponding 
transport skewnesses 7, Eq. ©, should be taken. Note 
that repulsive impurities generally lead to the opposite 
sign of 7. This can result in a partial suppression of spin 
Hall currents due to impurity compensation (where some 
acceptors are still present in an n-doped material and can 
localize electrons, leading to repulsive scatterers). 

Next we evaluate the contribution of the side jump [2]| 
to spin currents, where the wave function is laterally dis- 
placed during the scattering event. (This displacement 
does not modify the scattering angle measured at large 
distances, i.e., it does not affect the scattering cross sec- 
tion.) Side-jump currents were analyzed in detail for the 
anomalous Hall effect and we now relate the AHE to 
the SHE. In the AHE, a net polarization in combination 
with SO interaction at impurities, leads to electrical Hall 
currents, even in the absence of an external magnetic 
field. For the SHE, electrons are unpolarized in equilib- 
rium and we consider induced spin currents j M . Because 
A is small and spin relaxation is of order A 2 [25L |2(|, 
one can understand the e M -component of spin current as 
difference in particle currents of two spin species with 
polarizations ±e M . For non-interacting electrons, each of 
these species carries the anomalous Hall current of 
a system with density hah = \n and with spins fully 
aligned along the ±e M direction, and we can express the 
spin Hall current as 



J SH 
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(9) 



For the AHE, the side-jump contribution was found to 
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-2 n AH A (e 2 /ft) e M x E 9]. It results from SO 
corrections <5r to the velocity operator during impurity 
scattering. Nozieres and Lewiner clarified that this 
anomalous velocity <5r comprises two equal SO contribu- 
tions. The first is (5if = (i/K)[H, r] -v = (\/h)(VV xer) 
and becomes Sir = A(er x k) after the equation of mo- 
tion, k = —W/H, is taken into account. The second 



originates from the correction to the coordinate operator, 
the Yafet term <5rso = A(er x k) j2(|, and contributes as 
S2T = Sit. Note that 82Y leads to a factor of 2 which 
is often ignored. Heuristically, we can now understand 
the current J ah as follows. For scattering at an impurity 
with momentum transfer i5k, the lateral displacement it 
St = 2A(er x <5k). The anomalous Hall current is recov- 
ered from J A ^j = enSr/r, by using that the momentum 
dissipated per scattering event is Mk = — eEr. Using 
J A h and Eq. ©, we readily obtain the side jump 
contribution to the SHE. 



•fx 

jsj, 



= -2n\-e Kp - v E v 
n 



(10) 



The sum j 33 + j sj provides the total spin-Hall current. 

Mathematically, one could consider a model where the 
electron charge is cancelled by a uniform positive back- 
ground, with only small fluctuations in the potential V. 
In such a case, we could let t become arbitrarily large. 
Although the side-jump contribution to the spin Hall con- 
ductance is independent of t, the skew-scattering con- 
tribution, given by Eq. (JHJ, would grow with r. This 
growth would be cut off, however, when r becomes com- 
parable to the inverse of the SO splitting A due to the 
k 3 Dresselhaus term. If 7 is small compared to A /E-p, as 
would be the case if the potential fluctuations are suffi- 
ciently small, then the skew-scattering contribution will 
be smaller than the intrinsic contribution by a factor of 
order j Ep /A. 

We now evaluate the skewness 7 [Eq. ©] for a screened 
attractive Coulomb potential. For this, we make use of 
the long-established theory of single electron scattering 
by an atom J^]. We rescale the parameters to make 
a connection between the atomic Hamiltonian and the 
effective Hamiltonian [Eq. QJ] with V = —e~ qsT e 2 /er, 
effective mass m* , permittivity e, and screening length 
l/q s . We match V be setting the atomic number to 
Z = 1/e. Further, to match the SO interaction, we 
define an "effective" speed-of-light c* such that A = 
h 2 /4 (m*c*) 2 . Finally, when the sign of A differs from 
its vacuum value, we replace S(#) by —Sty) in Eq. © 
|27| . For GaAs, e = 12.4 and m* = 0.0665 m , thus 
c* w c/79 and a* Z w 1/21, with fine structure con- 
stant a* — e 2 /he*. When evaluating /($), screening of 
the long-range Coulomb interaction is required, other- 
wise the integrals in Eq. iJBJ diverge, e.g., r — > 0. For 
exponential screening and in second order Born approx- 
imation, J(#) is given by the Dalitz formula |22| (with 
the parameters used below and at large $, it agrees well 
with the exact solution for the unscreened potential). We 
assume Thomas-Fermi screening with inverse screening 
length q s = y // 3e 2 n/2eEp. The numerical evaluation be- 
low shows that the transport skewness 7 depends only 
weakly on the screening, implying that Eq. (JSJ is well- 
suited for evaluating the skew-scattering contribution. 
[Because of low sensitivity of 7 to screening, Eq. (JSJ is 
also in agreement with the common belief that the skew- 
scattering contribution is proportional to r.] When eval- 
uating £>($), we note that the main contribution comes 
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from large angles $ where the effects of screening are neg- 
ligible. We thus use for S(§ ) the exact expression for an 
unscreened potential @, |23 • Note that when evaluating 
the small parameter is (Za*) 2 = 4 |A| /(a* B ) 2 , with 
effective Bohr radius ag = h 2 e/m*e 2 , so we expect 1/500 
for the order of magnitude of 7. 

We now estimate 7 for GaAs and for the electron 
density n = 3 x 10 16 cm~ 3 reported in Ref. 01 i.e., 
Ef — 5.3 meV and qj 1 « 9nm. As a test, we first 
evaluate the longitudinal conductivity using an impurity 
density n,- = n and the Drude formula, a xx = e 2 nr/m* 
with t -1 = riiVp J dfl I (1?) (1 — cost?) and arrive at 
w 1.8 X 10 3 n _1 m -1 . This is within 10% of the 
experimentally observed conductivity at low voltages; a 
surprisingly good agreement, given that the rather small 
value kp£ ss 2 restricts the accuracy of the Boltzmann 
approach. Next we evaluate Eq. © and find 7 w 1/900. 
This value is rather stable and changes by less than 30% 
when q s or Ep increase or decrease by a factor of two. 

Next, we estimate the spin Hall currents. The 
measurements were performed at electrical fields E ~ 
20mV^m _1 where the conductivity increased to a xx ps 
3 x 10 3 Vl^tcT 1 due to electron heating. We assume 
that 7 is not very sensitive to these heating effects and 
we still use Eq. JSJ) but with the increased conductivity. 
For an electrical field E = xE x , we find both contri- 
butions to the spin Hall conductivity er SH = —jy/E x , 



namely ofg = -(j/2e)a xx ps 1.7 Cl~ 1 m~ 1 /\e\ ancl of? = 
2n\e/h ps — 0.8f2 _1 m _1 /|e|. In total, we arrive at the ex- 
trinsic spin Hall conductivity ofH ~ 0.9fi -1 m _1 /|e| pgj . 
The magnitude is within the error bars of the experimen- 
tal value of |of^, | w 0.5 f2 _1 m -1 /|e| found from spin 
accumulation near the free edges of the specimen 0, |29| . 

Bernevig and Zhang calculated intrinsic spin currents 
due to & 3 -Dresselhaus interaction ^(|. In contrast to 
the case of spin-orbit interactions linear in k, where ver- 
tex corrections lead to a cancellation of the spin current 
(l2l I30L l3lj. for the k 3 interaction this is not true; e.g., 
for s-wave scatterers the vertex corrections vanish iden- 
tically 32]. However, in the dirty limit tA/K ~ 10~ 2 , 
where A is the SO splitting A at the Fermi level and for 
experimental parameters of Ref. intrinsic spin currents 
were found to be small, only cr D « 0.02 fi -1 m -1 /|e| [l^ . 

In conclusion, we solved the kinetic equation including 
skew-scattering at impurities. We derived the contribu- 
tions to the extrinsic spin Hall effect and found reason- 
able agreement with experimental data. We acknowl- 
edge discussions with D.D. Awschalom, B.A. Bernevig, 
E. Dernier, Y.K. Kato, D. Loss, A.H. MacDonald, E.G. 
Mishchenko, A. Stern, R. Winkler, and S.-C. Zhang. This 
work was supported in part by NSF grants DMR-02- 
33773, PHY-01-17795, the Harvard Center for Nanoscale 
Systems, and DARPA. 
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Appendix A: ANGULAR INTEGRALS 

For completeness, we now show how the angular inte- 
grals in Eqs. (QJ and JSJ) can be evaluated. Without loss of 
generality, we choose the coordinate system k = (0, 0, k) 
and parameterize k' = k (cost/? sin sin ip sin cos &) 
since k' = k. We set a=(EC t +PxED t ) and Eq. |@Jl 
is obtained from 

J dft(k') / (#) (k - k') • a 

= J dfl I (#) [fc a z (1 — cosz9) — k a x cos ip sin $ 

— ka y sin <p sin #1 
= ka f dOI(tf) (1-costf), (Al) 

where we used in the last line that integration over dQ 
removes the terms in sin ip and cos ip and that ka z = 
k • a. Next, we note that n = k' x k/ |k' x k| = 
(sin <j), — cos <j>, 0) and insert into the left hand side of 
Eq. ©, 

f dO(k') / (i?) S (tf) (cr • n) (k • E — k' • E) C k 

= J dfi(k') I($)S(&) (tr* sin ^ - cr y cos ip) 

x [/c i? 2 (1 — cos d) — kE x cos <p sin $ 
— k E y sin sin d\ Ck 

= -~(ka x E y -ka y E x )C k J dQ(k') J(i?)S(0) sintf, 

(A2) 

using that also terms in sin y> cos p vanish after integra- 
tion and that J dp> cos 2 ip = f dtp sin 2 p = h J dip. From 
ka x E y — ku v E x = k • (cr x E), we recover the r.h.s of 
Eq. ©. 



